Linear conductance in Coulomb-blockade quantum dots in the presence of interactions 

and spin 



O 

o 

(N 
o 

Q 

m 



Y. Alhassid,^ T. Rupp,^ A. Kaminski,^ and L. I. Glazman^ 

^ Center for Theoretical Physics, Sloane Physics Laboratory, 

Yale University, New Haven, Connecticut 06520 
^Condensed Matter Theory Center, Department of Physics, 
University of Maryland, College Park, Maryland 20742-4111 
'^Theoretical Physics Institute, University of Minnesota, Minneapolis, MN 55455 

We discuss the calculation of the linear conductance through a Coulomb-blockade quantum dot 
in the presence of interactions beyond the charging energy. In the limit where the temperature 
is large compared with a typical tunneling width, we use a rate-equations approach to describe 
the transitions between the corresponding many-body states. We discuss both the elastic and 
rapid-thermalization limits, where the rate of inelastic scattering in the dot is either small or large 
compared with the elastic transition rate, respectively. In the elastic limit, we find several cases 
where a closed solution for the conductance is possible, including the case of a constant exchange 
interaction. In the rapid-thermalization limit, a closed solution is possible in the general case. We 
show that the corresponding expressions for the linear conductance simplify for a Hamiltonian that 
is invariant under spin rotations. 
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I. INTRODUCTION 

In the Coulomb-blockade regime of quantum dots, the 
conductance of the dot exhibits peaks as a function of 
the gate voltage, where each peak corresponds to the 
tunneling of one more electron into the dot. Of partic- 
ular interest are diffusive or ballistic chaotic dots, where 
the mesoscopic fluctuations of the conductance peaks and 
their spacings reflect the statistical nature of the groui. 
state wave function and energy of the isolated dot! 
The simplest model to describe Coulomb blockade is the 
constant-interaction (CI) model, in which the electrons 
occupy single-particle levels in the dot and the inter- 
action is described by an electrostatic Coulomb energy 
that is constant for a fixed number of electrons. In this 
model, the conductance near a Coulomb-blockade peak 
at temperatures that are large compared with a typical 
tunneling width can be derujed in a closed form using a 
master-equations approachBu This regime of sequential 
tunneling is relevant to most Coulomb-blockade experi- 
ments involving weakly-coupled dots. 

At sufhciently low temperatures, the conductance 
through the dot is dominated by elastic processes. How- 
ever, at finite temperature, the electrons in the dot may 
undergo inelastic scattering processes caused, for exam- 
ple, by electron-phonon or electron-electron interactions. 
In general, one may consider two limiting cases, depend- 
ing on the relative magnitude of a typical elastic tunnel- 
ing width Foi of an electron into the dot, and a typical in- 
elastic scattering width Fin of electrons in the dot. In the 
so-called elastic limit, Td 3> Tin, and inelastic scattering 
processes in the dot can be ignored. In the opposite limit, 
Foi Tin, inelastic scattering occurs on such a short time 
scale that electrons in the dot are effectively thermalized 
imniediately after an electron tunnels in or out of the 
dot.a This limit will be called the rapid-thermalization 



limit, which in the literature is also sometimes referred 
to as the inelastic limit. Statistical properties of the con- 
ductance peaks were studied withiii.the CI model in both 
elasticElu and rapid-thermaliKationDElfl limits, as well as 
for intermediate sit nations .113 

However, recent experiments in almost-isolated Cou- 
lomb-blockade dots clearly indicate that it is necessary to 
take into accoppt electron-electron interactions beyond 
the CI model.M For such interactions, the eigenstates 
of the dot are no longer characterized by Slater deter- 
minants. An important question that arises is how to 
calculate the linear conductance through such a dot. In 
this work we provide an answer in the limit of sequential 
tunneling. We make the general assumption that both 
the thermal energy and the typical excitation energies 
in the dot are small compared with the charging energy. 
This defines the regime of Coulomb blockade, where in 
the vicinity of the N + 1-st conductance peak only the 
manifolds of states with N and -I- 1 electrons in the 
dot contribute to the conductance (the manifolds with 
an electron number different from either N and A^ + 1 
are pushed away by the charging energy). Given the 
many-body eigenstates of the dot (with A^ and N + 1 
electrons), their energies, and the transition widths be- 
tween the states in the A^- and (A^ -I- l)-electron dot as 
an electron tunnels into the dot, we discuss the calcula- 
tion of the linear conductance. This is done by introduc- 
ing non-equilibrium probabilities of the dot to be in its 
various many-body states. These probabilities satisfy a 
set of rate equations describing the transitions between 
the many-body eigenstates. We consider stationary solu- 
tions in the linear-response approximation in both limits 
where inelastic scattering of electrons in the dot is weak 
(elastic limit) and strong (rapid-thermalization limit). In 
the elastic limit, the linearized rate equations cannot be 
solved in closed form in the general case, but there are 
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several important cases where a closed solution exits. 
Such cases include sufficiently low temperatures where 
only ground-state to ground-state, or ground-state to 
excited-states transitions (of N and -I- 1 electrons, re- 
spectively) contribute to the conductance. Another case 
where the rate equations can be solved in closed form at 
any temperature corresponds to a Hamiltonian for which 
the occupations of the orbital single-particle states are 
good quantum numbers. A particularly important such 
case corresponds to the universal HamiltonianElEj in the 
limit of infinite Thouless conductance. This Hamiltonian 
includes a constant exchange interaction in addition to a 
constant charging energy.EJ We also consider the rapid- 
thermalization limit, where we can always obtain a closed 
form expression for the conductance. 

Furthermore, when the Hamiltonian of the dot is spin- 
rotation invariant (as is the case of the Coulomb interac- 
tion), the expressions we find for the linear conductance 
can be simplified by carrying out explicitly the sum over 
the magnetic quantum numbers (in both the elastic and 
rapid-thermalization limits). 

The outline of this paper is as follows: In Section |l| we 
introduce the basic assumptions and notation. In Sec- 
tion HI we discuss the rate equations satisfied by the 



probabilities to find the dot in one of its many-body 
eigenstates with TV or iV -I- 1 electrons. We assume the 
elastic limit, where inelastic scattering of electrons in the 
dot is negligible. In Section IV we use linear-response 



theory to linearize these master equations in the station- 
ary limit. A general expression for the linear conductance 
in terms of the solution to the rate equations is derived 
in Section |^. We show that if detailed balance is sat- 
isfied for each pair of many-body eigenstates [in the iV- 
and {N + l)-electron dots, respectively], then the mas- 
ter equations and consequently the conductance can be 
solved in closed form. Furthermore, the corresponding 
expressions are shown to be simplified for Hamiltonians 
that are invariant with respect to spin rotations. Physi- 
cally relevant cases where such closed solutions exist are 
discussed in Section VI, and include the case of the uni- 



versal Hamiltonian in the limit of infinite Thouless con- 
ductance. In Section VII we solve the rate equations 



and derive a closed form for the conductance in the case 
where the inelastic scattering width is much larger than 
the elastic width (rapid-thermalization limit). Finally, 
in Section VIII we summarize the main results of this 
paper. 



to study the transport through the dot. Furthermore, 
we assume the Coulomb energy e^/C to be much larger 
than the thermal energy kT. In this limit, the linear 
conductance of the dot exhibits Coulomb-blockade peaks. 
When the gate voltage Vg is tuned in the vicinity of a 
Coulomb-blockade peak, the number of electrons in the 
dot can be either N or + 1. All manifolds with an 
electron number different from either iV or TV -f 1 are 
pushed away by the charging energy. 

The rate equations have been solved in the CI model 
for a single spin-degenerate leveH and for any number 
of single-particle levels.Q In the CI model, the two-body 
interaction is modeled by a constant charging energy of 
an object with capacitance C and charge Ne. Here, we 
derive and solve the rate equations in the presence of in- 
teractions beyond the charging energy, where the many- 
body states of the dots are no longer Slater determinants. 

We denote the many-body eigenstates of the dot with 
N electrons by i, and their respective energies by e^^'. 
Similarly, the states of the dot with TV + 1 electrons are 
labeled as j and the corresponding energies are e^'^'^^K 
Mn and M^j^i denote, respectively, the total number of 
the N- and (TV -|- l)-electron states considered. 

Of particular interest are cases where the dot's Hamil- 
tonian is spin-rotation invariant. This requires the ab- 
sence of a spin-orbit interaction such that the single- 
particle energies of the electrons in the dot are spin- 
degenerate. In addition, for spin-rotation invariance to 
be valid, the matrix elements of the two-body interaction 
may depend only on the orbital states but not on the spin 
indices, as it is the case for the Coulomb interaction. We 
then have 



[H, S^] = and [H, S^] = 



(1) 



The many-body eigenstates of the dot can be character- 
ized by the good quantum numbers of the spin S and its 
projection Sz = M. The remaining quantum numbers, 
in addition to the total number of electrons N, will be 
labeled by a. Thus the eigenstates of the dot with n 
electrons are given hy i = {a,S,M), and their respec- 
tive energies s^^J are independent of M. The states of 
the dot with -I- 1 electrons are similarly labeled by 
j = {a' , S' , M'), and their corresponding energies are 

^0^5' • example of a spin-rotation invariant Hamilto- 
nian is the recently proposed universal Hamiltonian.a't3 



II. MODEL 

We consider an almost-isolated quantum dot described 
by a Hamiltonian H that includes a one-body part and a 
two-body interaction. The dot is weakly coupled to leads, 
and we assume the limit kT ^ F, where F is a typical 
transition width of an electron from the leads into the 
dot. In this limit we can ignore the coherence between 
the dot and the leads and use a rate-equations approach 



III. RATE EQUATIONS 

We assume that a potential difference V is applied be- 
tween the two leads (source and drain) at temperature 
T. Fractions rji and rj^ of this potential difference, with 

VI + Vr- = 1 , (2) 

fall between the dot and the left lead, and the right lead 
and the dot, respectively. As a result a current / ffows 
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through the dot and various eigenstates of the dot can be 
occupied with different probabihties. We denote the non- 
equihbrium probabiUty for the dot to be in a particular 
n-electron state i by P^^^^ Since we consider situations 
in which the dot can only be occupied by or + 1 
electrons, we require the normalization 

^pW_,;^^(-+i)^l. (3) 

i 3 

We further assume that energy is conserved in the tun- 
neling between the dot and the leads, neglecting vir- 



tual transitions that are of higher order in the tunnel- 
ing widths. Denoting the transition widths between the 
iV-electron state i and the {N + l)-electron state j in- 
volving an electron tunneling into the dot from the left 
(right) lead by T\j (Tlj), we can write the following rate 

equations for the probabilities -Pj^^^ and pj^+^^: 



dP, 



dt 

{N+l) 

i 

dt 



E ^"^^'^[(1 - fli^3 + (1 - /&)ry - TSflA^ + fij^ij) fo^ au * , (4a) 

j 3 

E P^''\flA3 ^m,)- Pf^'^ E[(l - + (1 - /^)ry for aU J . (4b) 



The Fermi-Dirac function of an electron in the left (right) 
lead is evaluated at a suitable energy that can be found 
from energy conservation of the transition between states 
i and j, and is denoted by fj^ (fij)- In the presence of a 
gate voltage Vg, the energy of the n-electron dot in state 

i is given by e[^^ — Ne(Vg. Here C = Cg/C, with Cg 
being the dot capacitance with respect to the gate, and 
C being the full dot capacitance. Thus the energy of an 
electron (relative to the Fermi energy ep in the leads) in 
the left (right) lead that tunnels into an iV-electron dot 
in state i to form an {N + l)-electron dot in state j is 
given by 



time, 



(5) 



where ep = e(^Vg + ep is an effective Fermi energy. Taking 
into account the bias-potential drop between the dot and 
each lead, we have 

4 = /(£y +7?iey) and f,^ = f{e,, - VrcV) , (6) 

where the Fermi-Dirac function is defined by f{x) = [1 + 
exp(/3a;)]-i with P = 1/kT. 

We note that we can choose any values of rji and rj^- 
satisfying Eq. (^), but that the final result for the physical 
quantities (e.g. the conductance) must be independent of 
this choice (see the Appendix for a proof). 



IV. LINEAR RESPONSE 

To calculate the conductance, we are interested in find- 
ing stationary solutions of Eqs. (^, i.e. non-equilibrium 



probabilities P/^"* and P 



)(JV+1) 



that are independent in 



dp 



(N) 



dt 



= for all i . 



= for all j . 



(7a) 



(7b) 



Furthermore, to calculate the linear conductance, we are 
interested in calculating the current for small bias volt- 
ages V, where the equations can be solved in linear- 
response theory. To that end, we linearize Eqs. ^ in 
V. This involves an expansion of the Fermi-Dirac dis- 
tribution functions fjj, flj as well as the probabilities 

P^^^ and p^^^"*"^^ around their equilibrium values. The 
Fermi-Dirac functions are expanded up to first order as 

fl,^h,+meyn, and fl, ^ ~ V^eV fl^ , (8) 

where we used the short- hand notation /,j = /(sij). The 

probabilities Pj;^^ and pj^+^^ are expanded around their 

equilibrium values Pj^'^'^ and pj^^^-* which they would 
attain if the bias voltage would be zero and the dot would 
be in thermal equilibrium, 



P 



i 

(N+l) 



p(JV+l)r-^ 



(N+l) 



(9a) 
(9b) 



in terms of new variables 5'^^'' and '^^J^'^^K The equilib- 
rium probabilities can be expressed explicitly in terms of 
the eigenenergies e^^'^ and sf^^^^ of the dot with N and 
N + l electrons, 



P 



(AT) 



(10a) 



4 



and 



J 7. 



(10b) 



with the partition sum 



energy ep (at temperature T) in the leads, but because 
of the presence of charging energy only two values of the 
electron numbers are allowed. To find a set of equations 
for the new variables '^f^ and we substitute 

the expansions (||) and (^ into the r.h.s. of the rate 
equations (^ and equate the results to zero to ensure 
stationarity. Keeping only the terms linear in the bias 
voltage V, we obtain 



The equilibrium distribution is in principle grand- 
canonical with a chemical potential equal to the Fermi 



J 



3 

= Pl"^ ) f^^^3 + r-,) + E 4, (^ir^. - ^^r^,)] for aU ^, (12a) 

i 3 



(Jv+i)r«.T,(^+i) 



3 



Y.^~h3)(X\,+^lj)-Y.f^3^^3~^^3)\ iors^Wj. (12b) 



The above equations can be simplified utilizing the following three relations. From the definition ( p^ of the equilibrium 
probabilities, it follows that 



p{N+l) ^ p(Af)g-/3e., 



while the Fermi-Dirac distribution satisfies 



and 



-/3a: _ 



/(^) 



/'(x) = -/3/(x)/[l + exp(-/3:E)]. 
Using Eqs. (|l3|), (p^, and (|l5|), we arrive at the transformed rate equations 

E/^^[(rl,+r^,)(*f+'^-*W) + (r/irl^-r,,T^)] = O foralH, 

3 

E(i-/*^)[(rl,+r^,)(*f+'^-*f^) + (^irl,-^rr^,)] = o for all 



(13) 
(14) 
(15) 

(16a) 
(16b) 



Equations (^ represent a system of Mjsi -\- Mjv+i lin- 
ear equations for M^v -I- M^r+i unknowns. However, only 
Mtv -I- Afjv+i — 1 of these equations are linearly indepen- 



are de- 



dent. Hence, the solutions for ^^^^^ and 
termined only up to an additive constant. This constant 
is fixed by imposing the constraint 



E 



(17) 



r 



that follows from the normalization condition (^). In 
general, it is necessary to solve Eqs. ( [l^ ) numerically. 
These equations represent sets of detailed-balance equa- 
tions. However, it is possible to find analytic solutions 
in cases where a detailed-balance equation holds for each 
pair of states i and j individually, i.e. in cases where each 
term in the square brackets of Eqs. (f6|) vanishes. This 
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condition is equivalent to 



is given by the closed form 



(AT) 



{N+1) 



for any i and j . (18) 



Go 



e 

hkf 



(21a) 



We emphasize that there are M^Mn+i equations dl^ ) 
for only Mm + Mn+i unknowns, so in the general case 
it is not possible to satisfy (|l8|). However, there are a 
number of important cases where a solution does exist, 
as we discuss in Section VI. A solution to Eqs. (p^), when 
it exists, does not depend on temperature and leads to an 
expression for the conductance in which the dependence 
on temperature and on decay widths factorizes. 



V. LINEAR CONDUCTANCE 

The linear conductance G is defined as the ratio I jV 
in the limit ^ 0. The current through the dot can 
be expressed in terms of the quantities considered in the 
master-equation approach. In this framework, we can 
calculate for example the current / through the left tun- 
neling contact which equals the current through the dot 
and through the right contact. 



(19) 



or equivalently 



Gel = 



hkT ^ 



(21b) 



Notice that the parameters ryi^j. no longer appear in the 
conductance formula. The independence of the conduc- 
tance from the values of ryi^r, which holds true in gen- 
eral (cf . Appendix) , becomes apparent in this special case 
where an explicit formula can be given. For the conduc- 
tance in Eqs. (pl|), the temperature dependence of the 
contribution from each pair of states in the iV-electron 
and (iV -|- l)-electron dots is independent from the tran- 
sition widths. 

So far, the spin symmetries have not been taken into 
account in the derivation of the rate equations and the 
conductance. If this is done, the sum over the many- 
body states i and j can be further simplified. For two 
many-body states i = (a, S, M) and j = (a', S", M'), the 
associated partial decay widths V\a and FJ- are given by 



Substituting the expansions in Eqs. (||) and (^ into 
Eq. (p^), we find the elastic linear conductance 



Go 



(7V+1) 



(20) 



Although 771 appears explicitly in Eq. (pO[), the linear 
conductance is independent of 771. This must be so on 
physical grounds, and a formal proof is provided in the 
Appendix. 

In general, we have to solve Eqs. ( |T^ ) for "^f ^ and 
substitute their values in Eq. (^0|). However, 
in the special case where Eqs. (nsf) hold, the conductance 



v)':=T'^\{a'S'M'\^l{v,,)\aSM){' 



(22) 



where the operator ■!/'m(i"i) [V'm(i"r)] creates an electron 
with spin projection m at the left (right) point contact at 
ri (rr), and we have introduced overall coupling strengths 
Fg'"^ at the left and right point contact. The spin selection 
rules require m = M' — M and S' — [^±1/21 for the ma- 
trix element in Eq. ( p2| ) not to vanish. For a Hamiltonian 
which is invariant with respect to spin rotations, we can 
use the Wigner-Eckart theorem to factorize the matrix 
element in Eq. (^) into a Clebsch-Gordan (CG) coeffi- 
cient and a reduced matrix element that is independent 
of the spin-projection quantum numbers M, M', and m. 



{a'S'M'\^l{v)\aSM) = --j^={S M l/2m\S' M'){a' S'U\v)\\aS) 

= {-if-'' ( M' Vn M ) («'5'IIV'^(r)||a^) • (23) 

For completeness, we have given the relation in Eq. ( p3| ) both in terms of a CG coefficient (first line) and a Wigner 
3j-symbol (second line). For a given pair of many-body states i and j, the decay widths F^^- and V\j involve the same 
CG coefficient. In Eq. ([l8|), the CG coefficients in the ratio on the r.h.s. cancel and we are left with 

vpW _ ^ vA\{o^'S'W{vi)\\aS)\'' - ij,Tl\{a'S'U^{vr)\\aSr 

Tl\{a'S'U^{vi)\\aS)\^ + Tl\{a'S'U^{vr)\\aS)\^ ' ^ ' 
I 

Thus, if solutions for and exist that satisfy Eq. (|l^), they will be independent of the spin projec- 
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tion M or M', i.e. "^^^^ = ^^g' and ^1/^.^+^) = ^^^5^^^ Using the reduced widths in the expression (|l|) for the 
Defining reduced decay widths linear conductance, we have 

f:;5,o'5'=ro'''l(«'^'l!V'Uru)||a5)P, (25) 

conditions ( p^ ) can be rewritten in a form in which the 
magnetic quantum numbers disappear, 



aS.Q'S' aS,a'S' 

J 



aSa'S' MM'm aS.a'S' + aS.a'S' 



Taking advantage of the unitarity of the CG, J^Mm \{S M l/2m|S" Af')p — 1, we can write as our final result for the 



linear conductance 



^°^ = Wr ^ /(e^'S'.as)-^ —^^ • (28) 

aSa'S' ^ aS,a'S' + ^ aS,a'S' 



Instead of using the reduced widths defined in Eq. (25), 
it is possible to express the conductance in terms of the 
widths for the maximally-projected spin states (i.e. M — 
S and M' = S') 



aS.a'S' 



(29) 



The reduced matrix elements in Eq. ( p5D are non-zero 
only for S' = S'± 1/2 and, using the corresponding CG 
coefficients, we have 



rrs,a'5' = [2niax(5,5') + l] r'„^ 



S.a'S' 



(30) 



Ground-state transition at low temperatures 



At low enough temperatures, only the transition be- 
tween the ground states of the n- and {N + l)-electron 
dots [{N,S) {N + 1,5")] provides an important 
contribution to the conductance. It is then sufficient 
to consider the two ground-state manifolds {N, S) and 
(N + l, S') which are respectively {23+1)- and (2S" + 1)- 
fold degenerate. In this case, there is only; one equation 
( p^ ) and a solution can always be found ;Eil the two vari- 
are uniquely determined by 



ables ) and 



VI. EXPLICIT SOLUTIONS IN THE LIMIT OF 
ELASTIC SCATTERING 

In general, it is not possible to find a closed solution F'^ g, + g, 

for the conductance in the limit of elastic scattering, and 
Eqs. (Ilj) have to be solved numerically. However, in 
the following we list four important cases where explicit 

solutions to Eqs. ( p4| ) or more generally to Eqs. (|l^) exist and by the normalization condition (|l7|). The linear con- 
and the conductance is given by a closed expression. ductance then yields 



[2max(^,^0 + l]/(4^+^^-4^)-gF) T^s^s'^^S' 

{2S + 1) + (2^ + l)e-'3(4"^"-4"'-^F) Pg g, + T'g , 



where the widths Fg'^g, correspond to the maximally- [cf. Eqs. ( p9| ) and (^)]. According to Eq. (^2|), the func- 
projected spin states with M — S and M' — S' 
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tional form of the conductance (versus the effective Fermi 
energy) does not depend on tunnehng widths. In partic- 
ular, Gel is maximized when the effective Fermi energy 
ep = eCV^ + ep is tuned to 



_{N+1) 



-s 



kT ( 2S' + 1 \ 



B. Ground-state— to— excited-states transitions 

We consider only a single state in either the N- or 
the {N + l)-electron system and allow for any number 
of states in the other system. This is useful at low tem- 
peratures, when for an even number of electrons the low- 
est states in the spin sectors of S" = and S = 1 can 
be close in energy. The approximation of the previous 
paragraph may then be poor even at very low tempera- 
tures. For example, one might consider the transitions 
(iV, S = 1/2) ^ {N +1, a', S') with 5' = and 1. The 
variables associated to the {N + l)-electrons state are 
chosen as 



and the variable ^'L^'' is then fixed by the normalization 



77i.r 



S.a'S' 



(34) 



S,a'S' 



S,a'S' 



condition (|T^. Clearly, Eqs. (|3^) solve Eqs. (|24|) for all 
a' and the conductance has the form (|23). 



C. Systems with good orbital occupation numbers 

We define n\ — n\++n\^ to be the total occupation of 
the (doubly-degenerate) single-particle orbital A, i.e. the 
sum of the number of spin-up and spin-down particles in 
that orbital. Both n>,_|_ and n>,_ can take the value or 
1, hence n\ obtains the values 0, 1 or 2. We will show 
below that Eqs. (|2^) have a solution when all nx are good 
quantum numbers. For a given pair {a,S) and (a',S") 

only one term A = Aq in tp'^ir) = X)a '^^i^)'^\a ^^^^ then 
contribute to the respective matrix element 

{a'S'UHr)\WS) = ^VA(r)(a'5'|l4||a5) 

A 

= ijx„{r){a'S'\\alj\aS). (35) 

Clearly, if more than one orbital A contribute to the sum 
in Eq. (^), then the occupations cannot be good quan- 
tum numbers in the final state a' . Since the dependence 
of the reduced matrix element in Eq. (35) on the point 
contact at r is only through the wave function iJjxq (r) , wc 
have 



^ aS,a'S' 



(36) 



where 



r'/ = r|,'>,(r,,,)P 



(37) 



are the single-particle decay widths. 

We now argue that Eqs. (24) have a solution given by 



T\ + T\ 



nx{a,S) 



for TO = iV, iV + 1 , (38) 

where nx{a, S) are the single-particle level occupation 
numbers of the states with quantum numbers a and S 
[note that the label a includes all the good quantum num- 
bers {nx} = (ni,n2,...) plus any additional quantum 
numbers required to distinguish between states with the 
same occupations {nx} and spin S]. Since two manifolds 
of states {a,S) and (a',S") can only be connected by 
changing the occupation of a particular orbital Aq by a 
single particle, we have 



(JV) 
aS 



(JV+1) 
a'S' 



E 

A 



[nx{a, S) - nx{a , S')] 



Ao Ao 

Thus Eqs. (p6| ) can be indeed satisfied by 
conductance, we obtain the expression 



(39) 



For the 



G 



— y 

a'S' 



PisV(eaS.a'5') 



x\{a'S'\\al\\aS)\'- 



(40) 



where the orbital Aq depends on both the manifolds 
{a,S) and {a',S'). Of course, only manifolds that are 
connected by the addition of an electron to a single- 
particle orbital contribute to the sum in Eq. (|4C|); i.e. 
the occupations {nx{a' , S')} of the manifold (a', S") can 
only differ for one orbital and by one unit from the occu- 
pations {nxja, S)} of the manifold (q, S). An expression 
similar to (140) (but including explicit sums over the mag- 



netic quantum numbers) was discussed in Ref. 13 



An important example of a Hamiltonian in the class 
discussed here is the so-called universal Hamiltonian in 
the limit of the Thouless conductance gT ^ oo, 

^ = E ^>«L«A. + ^n^ - JsS' , (41) 



2C 



- aS,a'S' 



where S = X^a ^a with Sa = J2aa' 'Aa'^<^<^'^^<^' is the 
total-spin operator of the dot and <t is the vector of the 
three 2x2 Pauli matrices. Since [nx, Sx] = 0, the Hamil- 
tonian ( ^l|) is invariant under spin rotations and charac- 
terized by the good quantum numbers {nx},S, and M. 
The conductance through a weakly-coupled dot with a 
Hamiltonian ( pl| ) is then given by Eq. (pO|). 
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D. Systems with degenerate ground states 

In this subsection, we consider the case of low temper- 
atures, when only transitions between the ground states 
with TV and N + 1 electrons in the dot are allowed, and 
these ground states are M^r- and Mftr+i-fold degenerate, 
respectively. This case is more general than the similar 
case of Section VI A , since the degeneracy is not necessar- 



^'^^^^^ — "^''j"^^' +c, where c is a constant. The solution 
that satisfies the normalization condition (0) is the one 

with c = -(P(^) ^ + P^^^'^ Ej ^j)- While 

this constant c depends on temperature, it drops out in 
the final expression for the conductance, thus making the 
sum on the r.h.s. of Eq. (45) temperature- independent. 

The conductance in Eq. ( ^5| ) factorizes into two con- 
tributions. The sum over i and j is completely deter- 

1 r 

mined by the tunneling widths T^j and is independent 
of temperature and the effective Fermi energy (it is also 
independent of rji and yyi- as is shown in the Appendix). 
The prefactor of the sum on the r.h.s. of Eq. ( p5| ) does 
not depend on the tunneling widths T\'^ and contains the 
full dependence on the temperature and Fermi energy. 
This prefactor determines the functional dependence of 
the conductance peak on the gate voltage. Similar to 
Eq. (|3^) in Section VI A , the maximum of Gci is attained 
when the effective Fermi energy is tuned to 



ily the usual spin degeneracy, and the transition widths 

1 r • 

T^j can be arbitrary and are not necessarily related by 
a Wigner-Eckart theorem. Although we cannot write a 
closed solution for Gc\ in this more general case, the de- 
generacy of the levels participating in the transport pro- 
cess allows for significant simplification of the expression 
for Gel [as compared with the general expression given 
byEqs. (©, (|17|), and (|^]. 

The level degeneracy allows us to reduce the rate equa- 
tions (p^) to 



an, 



+ (r;irl, -77^^,)] -0 foralH, (42a) 



AN) 



2 Mn 



(46) 



j(rUry(* 



j 



+ (?;ir^,~77rrj,)]-0 for all J. (42b) 
As mentioned in Section IV, only + -/I/at+i — 1 of 



these Mat -I- M^r+i equations are linearly independent. 
The equation needed to determine a unique solution is 
given by Eq. (p7[). In the special case discussed here, 
the equilibrium probabilities for N and -I- 1 electrons 
are independent of the states i and j, respectively, and 
Eq. ( p7| ) simplifies to 



pW ^(N) ^ p{N+i) ^(.^+1) = 



j 



(43) 



where we have defined p(^) 

p{N+l 



p{N) 



pj^"*"^' for all j. In particular 
1 



M 



N ■ 



N+ie' 



-0{ei 



for all i and 



(44) 



where e^^^^ — e] for all i and e^^+i) = e 



[N+l) 



Using Eq. (^^ , the general expression (|^) for the 
ductance reduces to 



for all i . 



Here, ^I^ 



lie 



{N+l) 



£ 



{N) 



hkT Mn 

E 



X 



Afjv+ie 



-/3(e(« + i)-e(")-eF) 



(45) 



[N) 



and 'J' 



{N+l) 



are a solution to Eqs. (42) 



(N) 



These equations do not depend on the temperature, 
and we can therefore choose a private solution '5 
and 5'^^^^'' that is temperature-independent and deter- 
mined solely by the tunneling widths. The general so- 
lution of Eqs. (E2|) is given by ^ 



(N) 



M^'' + c and 



Eqs. (^2|) need to be solved n umeri cally. An exception 
is the case discussed in Section VI A , where the dot has 



spin-rotation symmetry and the degeneracy of the ground 
state corresponds to the various values of the spin pro- 
jection M. In this case Eqs. ( ^2|) can be solved in closed 
form [see Eq. (^l[)], and the linear conductance ( psf ) re- 
duces then to Eq. (||). 



VII. THE RAPID-THERMALIZATION LIMIT 

The above derivation of the linear conductance has 
been done under the assumption of dominantly elastic 
scattering, which is a good approximation at sufficiently 
low temperatures. Although the exact temperature de- 
pendence of inelastic scattering events is not well under- 
stood, they should become more relevant with increasing 
temperature, such that their width Fin will eventually be 
of comparable size to Fgi. We therefore consider in this 
Section the rapid-thermalization limit Fin S> Fd. 

We denote by P{N) the probability of the dot to be in 
a A^-electron state, and by P(j|A^) the conditional prob- 
ability of the dot to be in a particular many-body state 
i given the dot is occupied by N electrons. In general 



P-W = P{i\N)P{N) 



(47) 



In the rapid-thermalization limit, the conditional proba- 
bilities P(i|A^) and P{j\N + 1) are always given by their 
thermal equilibrium values P{i\N) and P{j\N + 1), re- 
spectively. It is then only the probabilities P{N) and 
P{N + 1) that obtain non-eqiulibrium values by a finite 
bias voltage V. These probabilities satisfy the following 
rate equations: 
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= P(Ar+i)^p(j|Ar+i)[(i_4)ri^. + (l_/^^.)ry-P(iV)^P(z|A^)[4 (48a) 

9t - p{N)Y,pm)[fi,Ti^+fi,Ti^]~p{N+i)Y,p{j\N+mi-n^^^^^ w 



dt 

dP{N+l 



As in the elastic case, we are interested in stationary 
solutions dP{N)/dt = dP{N+l)/dt = 0, where the rate 
at which electrons tunnel onto the dot (TV N + 1) is 
equal to the rate at which electrons tunnel off the dot 
{N + 1 ^ N). For a small voltage V, we can expand 
these probabilities around their equilibrium values P{N) 
and P{N + 1) to first order in V, 

P{N) ^ P{N) [1 + ey/?$(^)] , (49a) 
P(N + 1) ^ P{N + 1)[1 + ey/3$(^+i)] , (49b) 

where and are unknown variables. We ob- 



serve that while in the elastic limit there were Mjv + 
Mn+1 unknown variables, in the inelastic limit we are 
left with a considerably simpler situation of only two un- 
known variables. Consequently, an explicit expression 
for the conductance can always be given in the inelastic 
limit. 

In analogy to the elastic case, we proceed by substitut- 
ing linear expansions in Eqs. (^ and ( ^9| ) into Eqs. (|4^), 
and imposing stationarity. Collecting the terms that are 
independent of the bias voltage V, we obtain the usual 
detailed-balance equation at equilibrium. Keeping only 
terms linear in V, we arrive at 



PiN + Pm + i)[/3(i - + rj,)<i>(^+^) - 4,(^irJ, - vrTl^)] 

-P{N) F(^|iV)[/3/,,(^>^. + n^.)'i>^^^ + n,(jli^, - = . (50) 



Using the equilibrium relations P{N)P{i\N) ~ P^ ^^ and response current in Eq. (p^. We obtain 

PiN+l)Pij\N+l) = pI''+'^ 



together with Eqs. (13) 



(|TJ), and (|T^), we can rewrite Eq. ( pO| ) in the form 

V7^W/J($(-'+i)_cf>W)(rl,+rM 



-{r,iT\j - r,,Tl^)^ = 0. (51) 

The solution for <I)(^) and (up to an arbitrary 

additive constant) is then given by 

Sij Pi fiji^lj + ^ij) 

The additive constant is determined by the normalization 
condition P(A^)$(") + P{N + l)$(^+i) 0. 

To find the linear conductance Gthorm in the limit 
of rapid-thcrmalization, we substitute the expansions in 
Eqs. (H) and (B9|) in the general expression for the linear- 



Gthcon = ^ E ^ (a>(^) - $(^+1) - ,Orl, . (53) 



With the solution in Eq. (^ , we obtain a closed expres- 
sion for the rapid-thermalization conductance 



G 



therm ' 



Pi ^fij^\j^{j2kl Pk ''fkl^i 



hkT 



Srs P^ Vr-s(rJ,s + 



(54) 



As in the elastic case, we can exploit the spin symme- 
tries to perform explicitly the summation over the spin 
projections M and M', 

T.P'"'f^A= E ^"V(e.'S',.s)f 1;^..,., . (55) 

ij aSa'S' 

This yields as a final result for the conductance in the 
case of strong inelastic scattering 
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G 



the 



olSol' S' 



aSa' S' 



hkT 



aSa'S' 



(56) 



VIII. SUMMARY AND CONCLUSION 

We define 



(57) 



freedom of the dot and we showed the simphfications 
that occur when the dot's Hamiltonian is invariant un- 
der spin rotations. Both the hmits of dominantly elastic 
scattering and rapid thermalization are discussed. This 
work generahzes the solutions of Ref. ^ that were derived 
in the limit of non-interacting electrons (except for a con- 
stant charging energy). 



where Xij is a quantity that depends on the many-body 
states i and j of the N- and (iV-|- l)-electron dots, respec- 
tively. The equilibrium probabilities P^^^ and pj^+^-' to 
find the dot in states i and j are given by Eqs. ( |l0|) and 
(prj), while the energy difference is defined in Eq. (||). 

In the rapid-thermalization limit, the conductance of 
the almost-isolated dot is given by 



G 



therm 



(58) 



where F'^ {T\j ) are the partial transition widths between 
the states i and j involving the tunneling of an electron 
from the left (right) lead into the dot. 

In the elastic limit, a closed solution is not possible 
in the general case. The conductance can be calculated 



by solving the linear equations (|16[) for and 4" 



-AN+i) 



and then substituting in the expression (|20|). However, 
there are several important cases where a closed solution 
is possible, including the constant-exchange-interaction 
model (i.e. the universal Hamiltonian in the limit of infi- 
nite Thouless conductance). In these cases 
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Appendix 

Here we show that the linear conductance in Eq. ( |20| ) 
is independent of rj\ and T^r, the fractions of the bias po- 
tential difference between the dot and the corresponding 



leads. Defining the new variables 
can rewrite Eqs. ( p^ in the form 

E/'.[n + r^,)(*r^^-$f^) 



(N) 



(TV) 



ryi, we 



F"; 



)] = 



for all 



(60a) 



Gnl — 



hkT 



F' F"^- 
1. + Ff 



(59) 



When the Hamiltonian is invariant under spin rota- 
tions, these expressions can be simplified by carrying 
out explicitly the summation over the magnetic quantum 
numbers. The many-body levels of the N and {N + 1)- 
electron dots are now characterized by the quantum num- 
bers a,S,M and a',S',M' and the respective energies 
are independent of M and M'. Equations (^8|) and ( |59| ) 
are now valid with the reduced widths f J^g^,^, replacing 



aSMa'S'M' 



'S'M'' ^^'^ summation in Eq. (B 



the widths F 

carried over aS and a'S' only (but not over M and M'). 

In conclusion, we have solved the rate equations and 
found the linear conductance in the presence of interac- 
tions in the dot (beyond the charging energy). In par- 
ticular, we have taken into account the spin degrees of 



E 



(i-/..)[(rl,+F^)(vi/^ 



)-Fj^)]=0 
for aU j . (60b) 



The solution of Eqs. (^0|) is determined up to an additive 
constant. Since rji and ryr do not appear explicitly in 
Eqs. (pO[), we can choose a private solution $1^'' and 
^j^^^ is independent of rji and t],-. The general 

solution is then given by = $1^^ + c and = 

^j^^^^ -|- c, where c is a constant. In particular, the 
normalization condition ( pT|) can be satisfied by choosing 



(61) 
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The constant c depends on iji, but disappears in the final Expression ( |62| ) for the conductance shows clearly its in- 
expression for the conductance in Eq. (|2Cl|), dependence of 771 (and rj,-). 

Ge.= ^E^rV.($f^-*r^^)rl,. (62) 
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